% 创建时间：2025-06-21 23:57:39
% 作者：兪馳
clc;
clear;

A=[-3,-1,-1;
    0,-3,-1; 
    1,2,0]; 
j=sqrt(-1); 
A1=(expm(A*j)-expm(-A*j))/(2*j)
